Conserved Ordering Dynamics of Heisenberg Spins with Torque 
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We show that a torque induced by the local molecular field drives the zero-temperature ordering 
dynamics of a conserved Heisenberg magnet to a new fixed point, characterised by exponents z = 2 
and A ~ 5.15. Numerical solutions of the Langevin equation indicate that theories using a Gaussian 
closure are inconsistent even when the torque is absent. The torque is relevant even for quenches 
to T c , with exponents z — 4 — e/2 and A = d (where e = 6 — d). Indeed A is always equal to d for 
quenches to T c whenever the order parameter is conserved. 
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When a magnet is quenched from its disordered high 
temperature phase to its ordered configuration at low 
temperatures, the slow annealing of "defects" separat- 
ing competing domains, makes the dynamics very slow. 
The system organizes itself into a self similar spatial dis- 
tribution of domains characterised by a single diverging 
length scale which typically grows algebraically in time 
L(t) ~ t x l z . This spatial distribution of domains is re- 
flected in the scaling behaviour of the equal-time corre- 
lation function C(r, t) ~ f(r/L(t)). The autocorrelation 
function, C(0, ; 0, t) ~ L(t)~ x is a measure of the mem- 
ory of the initial configurations. The exponents z and A 
and the scaling function f(x) characterise the dynamical 
universality classes at the zero temperature fixed point 
(ZFP) 0. 

There has been a trend in recent years to compare the 
theories of phase ordering dynamics with numerical sim- 
ulations of Langevin equations. Comparison with experi- 
mental systems, such as magnets, binary fluids or binary 
alloys have to take into account the various 'real' features 
that might be relevant to its late time dynamics. For in- 
stance, theories of binary fluids have to include effects of 
hydrodynamics, while those of binary alloys have to in- 
corporate elastic and hydrodynamic effects. In the same 
vein, any comparison with the dynamics in real magnets 
has to include the effects of the torque induced by the 
local molecular field. 

In this article, we study the conserved phase ordering 
dynamics of a Heisenberg magnet in three dimensions in 
the presence of a torque (the corresponding nonconserved 
dynamics has been studied in Ref. Q). For quenches 
to T = 0, our Langevin simulation conclusively shows 
that the torque drives the dynamics to a new ZFP char- 
acterised by z — 2 and A ~ 5.15. This is confirmed 
by simple scaling arguments, which in addition gives the 
crossover time and the crossover exponent. Further we 
demonstrate that the approximate theories of conserved 
dynamics for vector order parameters based on the Gaus- 
sian closure scheme ||, are internally inconsistent even 
when the torque is absent, contrary to what has been 
assumed in the literature Q| . 

We investigate the dynamics following a quench to the 



critical point T c and show that the torque is relevant 
at the Wilson-Fisher fixed point. Using a diagrammatic 
perturbation theory, we show that z = 4 — e/2 at this 
new fixed point (where e = 6 — d) . We show to all orders 
in a perturbative expansion, that the autocorrelation ex- 
ponent A = d. This last result is true whenever the order 
parameter is conserved and can be understood from very 
general arguments [Q. 

The Heisenberg model in 3-dim is described by a vector 
order parameter <fi(r, t) and a free-energy functional, 
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This leads to the following dynamical equation (in 
dimensionless quantities), describing zero temperature 
quenches of the conserved Heisenberg model , 



dt \ 



(2) 



The first term on the right side is the dissipative force, 
while the second term is the torque generated by the local 
molecular field. The dimensionless parameter g ~ f2/F 
is the ratio of the precession frequency to the relaxation 
rate and is in the range 10 -3 — 10 for real magnets. 

We discretize Eq. @) on a simple cubic lattice (with 
size N ranging from 50 3 to 60 3 ) adopting an Euler scheme 
for the derivatives. We include the next-nearest neigh- 
bour contributions in calculating the second and higher 
order spatial derivatives which reduces lattice anisotropy 
and gives better numerical stability |5j. The space and 
time intervals have been chosen to be Ax = 2.5 and 
At = 0.2 which lead to stable results for the resulting 
coupled map. We have checked that slight variations 
of Ax and At do not change the results. Throughout 
our simulation we have used periodic boundary condi- 
tions. We prepare the system initially in the paramag- 
netic phase, where <j) is distributed uniformly with zero 
mean (angular brackets denote an average over this dis- 
tribution and over space). In our simulations we find 
that averaging over 5 initial uncorrelated configurations 
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gives us clean results (the statistical errors in C(r, t) are 
at most 3%). 

We compute the equal-time correlator C(r,t) = 
{(j){r,t) ■ cj)(0,t)}, the autocorrelation function, C(0,ti = 
0^2) = (0(0,0) • 0(0, t 2 )) and the energy density. Fig- 
ure 1 is a scaling plot of C(r,t) versus r/L(t) for various 
values of the parameter g, where L(t) is extracted from 
the first zero of C(r,t). Note that the scaling function 
for g = is very different from those for g > ; further 
the g > scaling functions do not seem to depend on the 
value of g. This suggests that the torque term drives the 
dynamics to a new ZFP. This is also revealed in the val- 
ues of the dynamical exponent z. In Fig. 2, a plot of L{t) 
versus t gives the expected value of z — 4 when g = 0. 
For g > 0, we see a distinct crossover from z — 4 when 
t < t c (g) to z = 2 when t > < c (g). The crossover time 
t c (g) decreases with increasing g. The same z exponent 
and crossover are obtained from the scaling behaviour of 
the energy density, e = ^ J dr ( ( V0(r, i) ) 2 ) ~ L(t)~ 2 . 



We provide a simple scaling argument to understand 
some of the results quoted above. On restoring appro- 
priate dimensions, the dynamical equation Eq.(0) can be 
rewritten as a continuity equation, d<f>(r,i)/dt — — V • j 
where the "spin current" 
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Using a dimensional analysis where we replace j a by 
the 'velocity' dL/dt, we find 



<tM 



dL = a_ 
dt L 3 ' " L 



(4) 



where M , a and T _1 are the equilibrium magnetisation, 
surface tension and spin mobility respectively. Beyond a 
crossover time given by t c (g) ~ (T/MqH) 2 ~ simple 
dimension counting shows that the dynamics crosses over 
from z = 4 to z — 2 in conformity with our numerical 
simulation. 




FIG. 1. Scaling plot of the equal time correlation func- 
tion for N — 50 3 . Observe the crossover in the scaling 
function f(x) as g changes from g = 0(o) to nonzero val- 
ues, g = 0.1(*), g = 0.3(A) and g = 0.5(D). 

The A exponent extracted from the autocorrelation 
function C(0, 0, t) also suggests a crossover to a new fixed 
point. The numerical determination of A is known to be 
subject to large errors (|], and so we have to go to very 
late times and hence large system sizes. For the 60 3 lat- 
tice we take data when 900 < t < 9000 and average over 
7 initial configurations. This time range is well within 
the scaling regime for C(r,t) when g = 0, 0.1, 0.2, 0.3. 
We compare the data for each g to A(t + t )~ x / z with A, 
to and A as fit parameters. The data show a crossover to 
a new ZFP beyond a time t c (g). We quote the values ob- 
tained for t > t c (g) — X(g = 0) w 2.2, X(g = 0.2) w 5.10, 
X(g = 0.3) w 5.17. 




FIG. 2. log-log plot ofi(t) versus t. At g = (o), z = 4 
(line of slope 0.25 drawn for comparison). At other values 
of g (g = 0.1(*), 0.3(D), 0.5(A)), z crosses over from 4 to 
2 (line of slope 0.5 drawn for comparison). Inset shows a 
scaling plot of y = L{t, g)jt x ^ versus x = tg 17 (see text), 
for g = 0.03(g), 0.05(*), 0.07(x), 0.09(D), 0.10(A). The 
scaling function asymptotes to a line of slope 0.25 (solid 
line) as x — > 00. 

The g dependence of the crossover time t c is also borne 
out by our numerics (Fig. 2(inset)) as we demonstrate 
below. Equation([|) suggests the scaling ansatz L(t,g) = 
t a s(tg^), valid for all g. This scaling form is governed by 
the 3 = ZFP, and so the scaling function s(x) should 
asymptote to a ^-independent constant So as x — > 0. This 
implies that a — 1/4. On the other hand, as x — > 00, 
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If the above proposal is true, then the data 



for L{t,g)lt x l A versus tg^ should collapse onto a scaling 
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curve for an appropriate value of <fi. We find that ^ « 1.7. 
Thus the crossover time, obtained when tg^ = 1, goes 
as t c ~ g^ 1 ' 7 , close to our earlier estimate. Note that 
our numerical estimate of the crossover exponent can be 
improved by introducing finite-time shift factors. 

Can the above results be understood within the class 
of approximate theories based on the gaussian closure ap- 
proximation |^Jl[]? The gaussian closure method consists 
of trading the order parameter <p(r, t) which is singular at 
defect sites, for an everywhere smooth field m(r,t), de- 
fined by a nonlinear transformation, <p(r, t) — a (m(r, t)) . 
Correlation functions are calculated making the single 
assumption that each component of m(r,i) is an inde- 
pendent gaussian field with zero mean at all times. The 
equal time correlation function takes the form H 
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where B(x, y) and 2*1(0, 6, c; z) are the Beta and hyper- 
geometric functions respectively and j(r, t) = (m(r + 
x,i) • m(x, i))/(m(x, t) 2 ). This scheme which works re- 
markably well for nonconserved systems || fails to give 
consistent results for our conserved dynamics, as we 
show, using a criterion developed by Yeung et. al. || 
for a conserved scalar order parameter. 
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FIG. 3. The spectral density 7O, t) at t = 3600 be- 
comes negative for < k/k m < 0.5 and for 1.5 < k/k m < 
3.0 (inset). 

We numerically evaluate the spectral density, the 
fourier transform of j(r, t). This is prone to numerical er- 
rors because of statistical errors in our computed C (r, t) . 
For instance, a numerical integration of J drC(r, t) pro- 
vides a nonzero value whereas it should be identically 
zero because of the conservation law. This is reflected 
in large errors in 7(k, t) at small k. We therefore adopt 
the following procedure. We fit a function Cf(x) to 



the equal time correlation function C(r, t) and use this 
to extract j(k,t) from the Eq(||). We observe (Fig. 3) 
that the spectral density, which should be a strictly pos- 
itive function of its arguments, becomes negative for 
k/k m < 0.5 ("f(k,t) is peaked at k m ) and in the range 
1.5 < kjk m < 3.0 ! This demonstration highlights the 
intrinsic inconsistency of the gaussian approach for con- 
served vector order parameters. 

We investigate the dynamics of the order parameter 
quenched to the critical point T c . At the Wilson-Fisher 
fixed point, u* = (8/ll)7r 2 e (e = 4 — d), the scaling di- 
mension of g is d/2 + 1 — z + 77/2, where z = 4 — 77 and 
r] = (5/242)e 2 implying that the torque g is relevant when 
d < 6. 

We calculate the dynamical exponents z and A at 
this new fixed point. This is done using a diagram- 
matic perturbation theory within the Martin-Siggia-Rose 
(MSR) formalism g. For our problem, the MSR gener- 
ating functional is, Z[h, h] = J X>(0)X>(0) exp (— J\<f>, 0] — 

H o[M + J °° dt J dk(h k ■ <?L k + h k ■ 0_ k ) where the MSR 
action is 



J[<l>,<l>] = J dt I dk^ i> k 



d, 



Tk 



,6F$ 



dkl ( ^-(k 2 - (k - k!) 2 )0 kl X k _ kl 



-Tk 2 c 



(6) 



Ho denotes initial distribution (gaussian with 
1 



and 

width Tq L and spatially uncorrelated) of the order pa 
rameter H = J dkf (4(0) • 0_ k (O)) @. 

Power counting reveals the presence of two different 
upper critical dimensions coming from the quatric term 
(d" = 4) and the cubic torque term {d 9 c = 6) in the action 
J. This implies we have to evaluate the fixed points and 
exponents in a double power series expansion in e = 4 — d 
ande = 6-rf §. 

The unperturbed correlation C k (<i,i2) and response 
C k (ti,t2) functions, and the bare u and g vertices are 
shown in Fig. 4. Again power counting shows that at 
d — 3, our perturbation expansion does not generate ad- 
ditional terms other than those already contained in J, 
i.e. the theory is renormalizable. However the pertur- 
bation theory gives rise to ultraviolet divergences which 
can be removed by adding counter-terms to the action. 

To remove these divergences, we introduce renormal- 
ization factors (superscripts R and B denote renor- 



malized and bare quantities respectively), 
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9 R = Z- X g 



L(0) = 
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It^B 
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Since the dynamics obeys detailed balance, the renor- 
malization factors Z and Z„ are the same as in statics. 
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Further the conservation of the order parameter forces 
ZZ = 1 to all orders. 

The new fixed point is given by the zeroes of the f3 
functions of the theory. The (3 functions, calculated from 
the Z factors, get contributions from all diagrams con- 



taining the primitively divergent diagrams T 

r (4) , T (Fig. 4). 



(2) 



r (3) 7 and 



The new fixed point, to one loop, is given by g* — 
±V192tt 3 £ + 0(£ 3 / 2 ), u* = (8/ll)7r 2 e + 0(e 2 ) (note u* 
does not change from its WF value to all loops) and the 
dynamical exponent z = 4 — e/2 + 0(e 2 ) [pi. 
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FIG. 4. Unperturbed (a) response function G k , (b) cor- 
relation function C£, and the (c) two bare vertices u and 

g. Wavy and straight lines represent the (j> k (t) and </>k(t) 
fields respectively, (d) Primitively divergent diagrams 
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The A exponent can be computed from the response 
function Gk(i, 0) = (</>k(0) • </>_k(i)) since this is equal to 
the autocorrelation function To _1 (<^k(i) 1 </'-k(0)), as can 
be seen from the first term in J on integrating by parts. 
The response function gets renormalized by 



G*(t,0) = Z- 1/2 G*(t,0). 



(7) 



The divergent contributions to Gb could come from 
two sources. Each term in the double perturbation se- 
ries could contain the primitively divergent subdiagrams 
p( 2 ) ; p(3) or y( 4 \ which we replace by their renor- 
malized counterparts. The other divergent contribu- 
tion could arise from the primitive divergences of the 
1-particle reducible vertex function r( 2 )(k, t,0), defined 
by G k (t,0) = /G k (i - t')rW(k,t',0)dt'. The superfi- 
cial divergence of the diagrams contributing to Gk(i,0) 
is D = V u {d - 4) + ^(d - 6) - 2 (where V u (V g ) is the 
number of u (g) vertices respectively). This is always 



negative for d < 6. This implies that G B (t,0) does not 
get renormalized and Z = 1. Consequently A stays at 
its mean-field value of d for this conserved Heisenberg 
dynamics both with and without the torque |]l0| . 

We have shown that the inclusion of a torque to the 
ordering dynamics of a conserved Heisenberg magnet, is 
relevant both for quenches to T = and T = T c . The 
new zero temperature fixed point is characterised by ex- 
ponents z = 2 and A w 5.15. We have provided scal- 
ing arguments to understand these exponents and the 
crossover. We have shown that the class of approximate 
theories based on the Gaussian closure scheme which 
had been constructed to understand this zero tempera- 
ture conserved dynamics, are inconsistent even when the 
torque is absent. On the other hand, the new critical 
fixed point is characterised by exponents z = 4 — e/2 and 
A = d (where e = 6 — d). Indeed A is always equal to 
d for quenches to the critical point whenever the order 
parameter is conserved. 
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